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Abstract 

We consider a singlet Dirac fermion with Peccei-Quinn(PQ) sym- 
metry as dark matter. A singlet complex scalar is introduced 
to mediate between dark matter and the SM through Higgs por- 
tal interaction and electroweak PQ anomalies. We show that a 
resonant annihilation of dark matter with axion mediation can 
explain the monochromatic photon line of the Fermi LAT data 
at 130 GeV by anomaly interactions while the annihilation cross 
section with Higgs portal interaction is p-wave suppressed. We 
discuss the interplay between the direct detection of the fermion 
dark matter and the collider search of Higgs-like scalars. We also 
present a ultra-violet completion of the dark matter model into 
the NMSSM with PQ symmetry. 



1 Introduction 



Dark matter is a dominant component of matter density in the universe, occupying about 
5 times larger than baryonic matter pp. It is known that dark matter(DM) interacts with 
the Standard Model particles gravitationally, while the property of dark matter such as 
mass and other interactions has not been known. Weakly Interacting Massive Particles 
(WIMP) has been a paradigm for cosmology, explaining the dark matter relic density by 
thermal freezeout with weak-scale mass and weak interactions for dark matter. Thus, 
it is expected that WIMP is detectable in underground direct detection experiments [2]. 
Recently, from the direct production of dark matter, the LHC has provided a new constraint 
on dark matter models, in particular, from the limit on spin-dependent cross section of 
dark matter [3]. 

While DM direct searches have imposed strong limits on the direct detection cross 
section, indirect searches from the cosmic gamma ray such as the Fermi Large Area 
Telescope (LAT) |4] has reached the sensitivity to DM annihilation cross section close 
to (av) = 3 x 10 _26 cm 3 s _1 , the one at DM freezeout, in the case of the purely s-wave 
annihilation [5j. The gamma ray constraints on the monochromatic photons coming from 
the DM annihilation have been reported from the Fermi LAT with 11 months and 24 
months data [6j[7]. A recent tentative analysis of the gamma ray data from the Fermi 
LAT [8],[9] has shown an indication for a gamma ray line at E 1 « 130 GeV at 4.6a local 
significance [9]. If the photon excess is explained by dark matter annihilating into a pho- 
ton pair, the observations hint at a dark matter mass of mx = 129.8 ± 2.4lJ 3 GeV and 
a partial annihilation cross section of (av)xx^-y-y = (1-27 ± 0. 32^28) x 10~ 27 cm 3 s _1 or 
(<jv)xx^- ry = (2.27 ± 0.57toil) x 10 _27 cm 3 s _1 , depending on the dark matter profile^. 
The results are very interesting, although it needs a further confirmation for excluding the 
possibility of instrumental errors or unknown astrophysical backgrounds. 

On the DM model building side, however, it is a nontrivial task to obtain such a 
large branching fraction into monochromatic photons from the annihilation of dark matter, 
because the relevant annihilation channels are loop-suppressed p~TVf2Tj . Furthermore, since 
the light SM degrees of freedom produced from the DM annihilation generates a continuum 
spectrum of gamma rays, one has to check if the continuum spectrum is consistent with the 
line spectrum. In this work, we propose a new dark matter model to explain the branching 
fraction into monochromatic photons as indicated by the Fermi LAT data [9] and discuss 
the gamma ray constraints in the model. 

We consider a fermion dark matter which carries PQ charge and couples to a complex 
scalar singlet. The complex scalar singlet does not couple to the SM directly. Instead, 
the real part of the complex scalar (CP-even scalar) mediates dark matter interactions to 
the SM Higgs boson by Higgs portal interaction, so does the imaginary part (the so called 
axion in our model) to the SM electroweak gauge bosons by PQ anomalies. We show that 

1 See also a different analysis of the Fermi LAT data [10] 
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the DM annihilation with axion mediation is a dominant channel in determining the relic 
density at freezeout by a resonance effect and it produces a monochromatic photon line 
by the DM annihilation into 77 or Z7 with sizable branching fractions at present. On 
the other hand, the DM annihilation with CP-even scalar mediation can contribute to the 
total cross section at freezeout comparably to the one with axion mediation while becoming 
p-wave suppressed at present. As a consequence, we show that there is a parameter space 
that explains the observed Fermi gamma ray line at E~ ~ 130 GeV and satisfies all the 
phenomenological constraints such as gamma ray constraints, DM direct detection and the 
Higgs-like scalar search at the LHC, etc. 

Although the introduction of extra DM annihilation coming from the CP-even singlet 
scalar reduces the gamma-ray line to the central values as given in Ref. [9], extra con- 
tribution does not have to be sizable, given the systematic and statistical errors of the 
DM annihilation cross section reported. In other words, the predicated gamma-ray line of 
our model without CP-even scalar mediation is still consistent with DM annihilation into 
two photons (at about 2 sigma level with the NFW dark matter profile [9]), so is it with 
DM annihilation into Z7 for any dark matter profile. But, we have included the CP-even 
singlet sector in our DM discussion as it naturally appears in the PQ completion of the 
axion-like scalar. 

The paper is organized as follows. First, we begin with a description of our model with 
PQ symmetric dark sector. Then, we present the results of the DM relic density and the 
production mechanism of gamma ray line and discuss the direct detection constraint on 
the model. We proceed to present a ultra-violet (UV) complete model for generating the 
electroweak anomalies by the electroweak PQ axion in the context of the NMSSM. Finally, 
a conclusion is drawn. There are two appendices providing the details of the decay rates of 
the CP-even scalars and the axion, and the cross sections for dominant DM annihilation 
channels with axion and CP-even scalar mediations. 



2 The Model 

We consider a Dirac fermion dark matter x that is charged under U(1)pq global symmetry. 
When a complex singlet scalar S couples to the fermion dark matter, it mediates the dark 
matter interactions to the Higgs boson H by the Higgs portal interaction and the Standard 
Model electroweak gauge bosons A 1 ^ by anomalies. The SM Higgs doublet does not have a 
PQ charge because of the Yukawa couplings [22] ■ We define PQ transformations on dark 
matter and the mediator as x ~ ► e l75Q x and S — > e~ 2ta S, respectively. 
We specify the relevant action for dark matter to be the following, 

C = ixi^x + \d,S\ 2 + \D^H\ 2 - V(H, S) - \fLF*» + Ant (1) 
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with 



V(H, S) = \ H \H\* + \ S \S\* + 2X HS \S\ 2 \H\ 2 + m 2 H \H\ 2 + m 2 s \S\ 2 - (±m' 2 S 2 + h.c.)(2) 
C mt = \ x (SxPlX + S*xPrX)+ E ^aF; u F^ (3) 

1=1,2,3 nVs 

where F^ v = |e M y p(J F pcr , v s = \/2 (5*) is the axion decay constant and we need |A X | < 0(1) 
for a valid effective theory for dark matter with mass M x = X x v s /\^2, and the constant 
parameters Cj depend on the anomalies generated by the axion-like scalar. If dark matter 
has a fixed mass as hinted by the tentative result from the Fermi LAT data, the dark 
matter coupling A x is determined by the axion decay constant v s . We take v s to be larger 
than the axion mass such that the anomaly loops with heavy fermions having masses of 
order v s are well approximated by the dimension-5 interactions between the axion and the 
SM electroweak gauge bosons. 

The interactions of a complex scalar field S = (s + ia)/ a/2 containing a as the imaginary 
part generates the DM mass and DM-axion coupling in the above action. When the 
anomalies with nonzero q are generated by the couplings of the complex scalar 5* to extra 
heavy fermions with SM charges, the axion mediates the dark matter interactions to the 
SM gauge bosons only through anomalies^). We assume that our axion does not couple to 
colored fermions so it can be called the electroweak axion. Otherwise, the DM annihilation 
into a gluon pair would be too large to give a sizable branching fraction into photons. Since 
the anomaly interactions are model dependent, we treat them to be arbitrary parameters. 
We postpone the discussion on microscopic models to a later section. 

For the PQ mechanism to work for solving the strong CP problem, the QCD anomalies 
must be generated by the invisible axion that couples to extra heavy quarks [22J. We 
note that the axion-like scalar a gets a PQ-breaking mass m' s from a higher dimensional 
interaction with a PQ breaking scalar $ containing the invisible axion with ($) = F a : for 
V = -2MT$ A S 2 + h.c, we get m' s = jfe. Thus, for F a ~ 10 10 GeV, which is within the 

invisible axion window, 10 9 GeV < F a < 10 12 GeV, the Planck suppressed term generates 
a weak-scale mass for the axion-like scalar. Furthermore, the invisible axion can constitute 
the dark matter relic density too, but we assume that it is subdominant. On the other 
hand, if the soft PQ-breaking mass violates CP, a mixing between CP-even and -odd scalars 
could lead to too large branching fraction of the DM annihilation into the SM particles. 
In order for the axion to couple to the SM only by anomalies, we assume that the induced 
PQ breaking mass does not violate CP. 

2 We note, however, that in the SM with two Higgs doublets H u and Hd, the Higgs doublets can 
have nonzero PQ charges [53] such that a renormalizable coupling, SH u Hd or S 2 H u Hd, could be written 
depending on the PQ charge assignments. In this case, the axion could mix with the pseudo-scalar Higgs, 
acquiring a tree-level coupling to the SM particles. As will be discussed in a later section, such couplings 
can be chosen to zero in a PQ supersymmetric NMSSM. 
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Figure 1: Feynman diagrams for dark matter annihilation with axion and CP-even scalar. 

After minimizing the potential in eq.((2]), the VEVs of the singlet and the Higgs doublet, 
v s and v, are determined as [21] 

2 ^Hsm 2 H - \ H {m 2 s - mf) 
2 A Hg (m| - ml) - X s m 2 H 

The conditions for a local minimum are \ijs m2 H ~ — m'g) > 0, Xnsi'^'s ~ m 's) ~ 

^s m H > and \ S \ H — A^ 5 > 0. Expanding the scalar fields around the vacuum as 
S = (v s + s + id) j y/2 and H T = (0,v + h)/ a/2 in unitary gauge, the obtained mass matrix 
for CP-even scalars can be diagonalized by the field rotation, 

s = cos 6 s + sin 6 h, h = — sin 6 s + cos 9 h (6) 

with 

tan 26 = - — — -, (7) 

X H v 2 - X s vi 

and the mass eigenvalues are 



ml 2 = \ H v 2 + X s v 2 s =f J(X s v 2 - X H v 2 f + AX 2 HS v 2 v 2 



Due to the singlet-Higgs mixingfl, the real scalar also mediates the dark matter interactions 
to the SM particles as the Higgs boson does. The axion mass is just given by the PQ 
breaking mass as m a = m' s . In our model, we can trade four independent Lagrangian 
parameters in the CP-even scalar sector to four physical parameters by eqs. (4)-(8): v s , m 12 
and the mixing angle 9. 



3 Recently, even a small singlet-Higgs mixing can lead to a sizable threshold correction to the Higgs 
quartic coupling at tree-level, solving the vacuum instability problem of the Higgs mass lighter than 
130 GeV in the SM [25]. 
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Figure 2: Parameter space of m a vs m 2 consistent with WMAP 3a band of the relic density 
(Blue). Left: Parameter space of the present partial annihilation cross section into a photon 
pair within 4 — 8% (Red) and larger than 8% (Green) and c\ = c 2 = 1. Right: Parameter 
space of the present partial annihilation cross section into Z'j within 8 — 13% (Red) and 
c\ = 0.2, C2 — 1. The mass of the Higgs-like scalar is chosen to ni\ = 125 GeV. 

3 Dark Matter Constraint and Fermi Gamma Ray 
Line 



In this section, we consider the constraint coming from the DM relic density in our model 
and show how a monochromatic photon line is produced by DM annihilation as observed 
by the Fermi LAT. We also discuss the interplay between DM direct detection and Higgs 
search at the LHC in the interesting region that explains the Fermi gamma ray line. 



3.1 Dark matter annihilation cross section 

Dark matter annihilates into an SM pair through both the axion and the real scalar partner 
of the axion. There are four channels with s-channel axion in our model as shown in Fig. 1: 
XX ~ >* a ~ >* 77> ZZ,W + W~ . All the channels are s-wave so their annihilation cross 
sections are little changed since the freezeout. The first two channels lead to monochromatic 
gamma lines from dark matter annihilation. If the branching fraction of the first two 
channels is sizable and DM annihilation occurs due to the resonance at m a ~ 2M X , it is 
possible to explain the tentative gamma ray line observed at Fermi LAT with the cross 
section required to explain the dark matter relic density. We denote the annihilation cross 
section coming from axion mediation by (av) a . 

On the other hand, dark matter can also annihilate by the mixing between the CP-even 
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singlet scalar and the Higgs boson. If kinematically allowed, there are four channels with 
s-channel CP-even scalar annihilating into the SM particles as shown in Fig. 1: XX ~^ 
s — > f f, ZZ, W + W~, hh. As shown in appendix B, all the channels are p-wave suppressed. 
So, even if they can be relevant for generating the correct relic density at freezeoutEl, they 
become suppressed at later time. We denote the annihilation cross section coming from 
CP-even scalar mediation by (crv) s . Moreover, singlet self-interaction and/or Higgs portal 
interaction could lead to extra annihilation channels at tree level: XX ~* ss > s h, hh, aa, as 
in the interaction basis. All the additional channels except XX ~~ >* as are p-wave suppressed. 
In particular, the first three channels into CP-even scalars can be comparable to the loop- 
induced channels with axion mediation at freezeout while they do not affect the branching 
fraction of the annihilation cross section into monochromatic photons. On the other hand, 
the XX ~^ aa channel is kinematically forbidden even at freezeout, close to m a = 2M X 
that we are interested in for gamma ray line. The XX ~~ >* as channel is s-wave so the 
cascade annihilation of the CP-even scalar into the SM particles could produce too many 
intense secondary photons [5] . So, we forbid the XX ~~ >* as channel kinematically by taking 
m a + mi j2 > 2M X and search the parameter space that is compatible with the relic density. 
Therefore, the total annihilation cross section is given approximately by the addition of 
the axion- mediated and scalar-mediated cross sections at freezeout, 



Then, from the velocity times cross section of dark matter annihilation, av = a + bv 2 , 
the dark matter relic density is given by 



where the freeze-out temperature gives x F = M x jTp « 20 and g* s (xp) is the number of 
the effective relativistic degrees of freedom entering in the entropy density. The DM relic 
density measured by WMAP is Q DM h 2 = 0.1123 ± 0.0035 at 68% CL pQ. 

In Fig. [2J on the m a — m 2 mass plane, for the fixed DM mass to M x = 130 GeV or 
145 GeV, we show the parameter space satisfying the WMAP bound on the relic density 
at 3a and the branching fraction of the DM annihilation into a photon pair within 8% or 
into Z^i within 8 — 13% at present. We have used the Feynrules [27] and micrOmegas [28] 
for the numerical analysis. There are four separate lines satisfying the relic density. Each 
pair of lines parallel to the «2 axis gets close to the axion resonance due to a small decay 
width of the axion. The necessary tuning for the axion resonance to produce the correct 
relic density within the WMAP 3a band is ~ 0.008 for m a = m° a + Sm a , that is, 
5m a ~ 2 GeV for m° a = 2m x = 260 GeV. Each pair of lines parallel to the m a axis are 
separated wider from each other due to a larger decay width of the CP-even scalars. We 

4 Dirac fermion dark matter with a CP-even singlet scalar mediation was previously considered |26j . 



(<7t>)|f r ~ (crf)a|fr+ (<™)s|fr- 



(9) 



n v h 2 



2.09 x 10 s GeV" 



(10) 



Mpi y/g* s {x F ) (a/x F + 3b/x 2 F ) ' 
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note that there is a parameter space in red in both figures where the gamma ray line at 
130 GeV can be obtained with the observed intensity as will be discussed in more detail 
later. 



3.2 Gamma ray line spectra and differential flux of photons 

The spectrum of the 77 line is a delta function at M x . On the other hand, the photon 
spectrum per annihilation with the process XX ~~ >* %1 depends on the mass and width of 
Z-boson as follows [16] . 

d_Nl _ AM X M Z T Z 
dE 1 /i/ 2 
where V z is the decay width of the Z-boson and 

, (M 7 \ , /4M^ - Ml\ 

h s ten -.(_i) +taa -i(_ fe7 »), (12) 

f 2 = (4Ml-4M x E.,-Mlf + T\Ml. (13) 

The differential photon flux produced by dark matter annihilation and integrated over 
the region of angular size AQ is computed as [16] 

dE^ ] ~ ^uqdR; JAn (14) 



with 



^ = 2w(av) r ,S(E,-M x ) + (av) Zy d ^., (15) 
) = J_/ JW , JW= / (16) 

JAQ Jlos r V P© 7 

Here, we note that p(x), p© = 0.3GeV/cm 3 and r Q = 8.5 kpc denote the dark matter 
density at a location x with respect to the GC, its value at the solar system and the 
distance between the Sun and the GC. And the coordinate s spans along the line of sight, 
making an angle ip with the direction of the GC. In order to compare to the spectrum 
observed by Fermi LAT, we should take into account a 10% energy resolution so we can 
replace the actual spectrum by a Gaussian distribution function centered at the photon 
energy. 

There are models for the dark matter density distribution in our galaxy. The Navarro 
Frenk and White (NFW) profile is often used for indirect searches, 

PNFwir) = Pa — -2 (17) 
f + f 
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Figure 3: Branching fraction of the partial cross sections for DM annihilations with axion 
mediation vs Ci/c 2 . Dark matter mass is fixed to give the dominant photon line. 



On very small scales, the optimal fit to simulated DM halos is provided by the Einasto 
profile, 

p(r) = p exp[-^((£f-l)] (18) 

where a = 0.17 and R = 20kpc. The Einasto profile is shallower than NFW at very small 
radii. 



3.3 Monochromatic gamma ray from fermion dark matter 

In this subsection, we consider the annihilation of the fermion dark matter with axion 
mediation and show that the monochromatic photons with sizable branching fraction can 
be produced. 

Suppose that the gamma ray line at E 1 = 130 GeV hinted by the recent analysis |9J is 
explained by dark matter with M x = 130 GeV annihilating into two photons. Then, there 

is another peak at £7 7 = M x ll — — 114 GeV. On the other hand, if the observed 

gamma ray line is explained by the DM annihilation into one photon, we need to choose the 
dark matter mass to M x = 145 GeV and another peak coming from the DM annihilation 
into two photons is at E 1 = 145 GeV. 

The ratio of the annihilation cross sections into Z7 and 77 in our model is determined 

by 

_ {(tv) Zi _ (c 2 a 2 - Ciai) 2 sin 2 (2gH/) / Mf y 
r ~ (av)^ ~ 2( Cl cti cos 2 6 W + c 2 a 2 sin 2 9 W ) 2 \ 4M 2 / ' [ ' 
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Figure 4: Partial annihilation cross sections vs ratio of CP-even scalar to axion-mediated 
cross sections at freezeout. Dark matter mass is fixed to give the dominant photon line. 
Left: WW, ZZ, 77, Z7 from top to bottom. Right: WW, ZZ, Z7, 77 from top to bottom. 



For ci = C2, which is the case where two Higgsinos generate electroweak anomalies in the 
NMSSM as will be discussed in the later section, the above ratio becomes r ~ 0.54 for 
M x = 130 GeV. Then, the intensity of the one-photon line coming from Zj is then 0.27 
times that of the 2-photon line. Although the one-photon peak may be resolved from 
the two-photon peak separated by 16 GeV, the suppressed one-photon line may not be 
significant due the background while the 2-photon line can explain the observed Fermi 
gamma ray lineH. On the other hand, for c\ = 0.2c2, the cross section ratio becomes 
r ~ 2.9 for M x = 145 GeV. Then, the intensity of the two-photon line coming from 77 
is 0.69 times that of the 1-photon line. But, in this case, the 1-photon line fits the Fermi 
gamma-ray line worse than the case with 2-photon line dominance [29J. The branching 
fraction for the annihilation cross sections is depicted in Fig. [3] for M x = 130 GeV and 
M x = 145 GeV, depending on ci/c 2 . 

In order to explain the observed Fermi gamma ray line by dark matter annihilation, one 
needs the cross section of dark matter annihilating into a pair of monochromatic photons 
to be (crv), n = (1.27 ± 0.32±g;i|) x lO^cmV 1 for the Einasto profile and (aw) 77 = 
(2.27 ± 0.57±g;|?) x lO^cmV 1 for the NFW profile, that is, Br(xx ^ 77) ^ 4 - 8% 
for thermal dark matter |9J. In our model, the DM annihilation with axion mediation 
is s-wave so it alone would have produced a more intense gamma ray line than observed 
because the branching fraction into photons is rather large as shown in Fig. [3J But, the DM 

5 It has been shown that two photon lines are shown to fit the gamma-ray data equally well [25]. In 
particular, if (av) z-y / \2{av) 1 ~ l ) < 1, the region around dark matter mass 130 GeV remains the best fit. 
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annihilation channels with CP-even scalar mediation can give a sizable contribution to the 
total cross section at freezeout while they do not affect the branching fraction into photons 
at present because of a p-wave suppression. Thus, as (av) s = at zero temperature, the 
total annihilation cross section at present is less than the one at freezeout so we can reduce 
the partial annihilation cross section into photons at present. We can write the partial 
cross sections at present in terms of the total annihilation cross section at freezeout, 



We plot the present partial annihilation cross sections with respect to the ratio of scalar- 
mediated to axion- mediated cross sections at freezeout in Fig. HI For c\ = C2, the DM 
annihilation into two photons is dominant to explain the Fermi LAT peak and we get 
Br(xx — > 77) ~ 0.14. In this case, we need the annihilation cross section with CP-even 
scalar mediation at freezeout to be (<rv) s / (crv) a = 0.8 — 2.3 for (crt>) 77 /(10 _27 cm 3 s _1 ) = 
1.27 — 2.27. On the other hand, for c\ = 0.2c2, the DM annihilation into one photon can 
explain the Fermi LAT peak and we get By(xx ~~ %l) — 0.13. So, we need a smaller 
annihilation cross section with CP-even scalar mediation as compared to the case with two 
photons: (crv) s /(av) a = - 0.54 for (^^/(lO^cmV 1 ) = 2.54 - 3.9. 

We note that the partial cross section into WW is the largest, being about 5 times 
larger than the one for the dominant photon channel. But, it satisfies the current limit 
of about 10~ 25 cm 3 s _1 , coming from the gamma ray emission of dwarf spheroidal galaxies 
observed by the Fermi LAT [5]. Furthermore, PAMELA has measured the spectrum of 
cosmic anti-proton flux below 180 GeV, which is consistent with the background [30] • In 
our model, the anti-proton flux coming from the hadronic decays of WW, ZZ with the 
annihilation cross sections of about 60% and 20%, respectively, is well below the measured 
one by one or two orders of magnitude [31] . 

4 Direct detection and LHC Higgs search 

As discussed in the previous section, although the DM annihilation with axion mediation 
gives a sizable branching fraction of the annihilation cross section into monochromatic 
photon(s), the level of the gamma ray line excess of the Fermi LAT data requires a sizable 
contribution of the CP-even scalar mediation to the total annihilation cross section at 
freezeout. This is achieved by a sizable mixing between the CP-even scalar singlet and 
the Higgs boson. Then, the same mixing parameter determines the direct detection cross 
section of dark matter. 

From the scattering process with the fermion dark matter, the spin-independent cross 
section for dark matter with nuclei is at tree level given by 



((Tv)xY 



Br(xx -> XY) 



•M fr , Bi(xx^XY) 



(<7V)XY 

(crv) a 



(20) 





(21) 
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Figure 5: Parameter space of the mixing angle of CP-even scalars vs the mass of singlet- 
like CP-even scalar. The region consistent with WMAP 3a band of the relic density is 
bounded by blue dotted lines while the blue solid line is the central value within 3a. Left: 
Parameter space of the present partial annihilation cross section into a photon pair within 
4 — 8% (Red region bounded by red solid lines) and less than 4% (Pink) and c\ — c-i — 1. 
Right: Parameter space of the present partial annihilation cross section into within 
8 — 13% (Red region bounded by red solid lines) and c\ = 0.2, c 2 = 1. The mass of the 
Higgs-like scalar is chosen to nil = 125 GeV. In both figures, gray region is excluded by 
XENON 100T and dark gray region is disfavored by the electroweak precision data. 

where m r = m N M x / {m N + M x ) is the reduced mass, m N is the nucleon mass, f^, 
parametrizing the Higgs-nucleon coupling, is given by the sum of the light quarks (/l) 
and heavy quarks (/#) as /jy = Yl II + 3 x -^fn [22], and m 12 are physical CP-even scalar 
masses given in eq. OH]). For instance, the direct detection bound is cr^_ N < 10~ 8 pb for 
M x = 130 GeV [2]. We note that the spin-dependent cross section with axion exchange is 
velocity-suppressed so the bounds coming from IceCube [33] and Super-Kamiokande [31] 
does not constrain our model. 

On the left plot of Fig. [5] in the parameter space of the mixing angle and the mass 
of the singlet-like CP-even scalar, for M x = 130 GeV, we considered the WMAP 3cr band 
for the relic density and the branching fraction of the annihilation cross section into a 
photon pair less than 8% and also show the limits from DM direct detection as well as 
the electroweak precision data (EWPD) at 95% C.L. [2HI35]. We have also shown the 
case with M x = 145 GeV on the right plot of Fig. [5] where both the WMAP bound and 
the gamma-ray line are obtained in all the parameter space away from the CP-even scalar 
resonance at m 2 ~ 2M X and the direct detection bound is little changed as compared to 
the case M x = 130 GeV. We can see that there is a parameter space with the mixing angle 
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smaller than about 0.4 — 0.5 that explains the Fermi gamma ray line at E 1 = 130 GeV 
and is compatible with all the phenomenological constraints. As shown in both plots of 
Fig. El the region below m 2 ~ 180 GeV is not constrained by direct detection because 
there is a cancellation effect in the scattering cross section due to the opposite signs in the 
amplitudes with CP-even scalars. In the range of the singlet-like CP-even scalar masses 
away from the cancellation zone, the mixing angle is constrained to be smaller than about 
9 = 0.4 by DM direct detection. We note that the EWPD does not give a stronger bound 
that the direct detection, in the region with the correct relic density and the Fermi gamma 
ray line. But, the former gives a stronger bound than the latter for m 2 < 260 GeV and 
m 2 > 350 GeV. 

We note that if the singlet-like CP-even scalar gets mass outside the range, 122 GeV < 
m 2 < 128 GeV, and there is no additional decay mode other than the ones of the SM 
Higgs boson, the LEP and LHC limits on the Higgs-like couplings will apply. The LEP 

restricts the Higgs-like coupling to Z-boson, £ 2 = (ghzz/ 9hzz^ to be less than 0.5 below 
m 2 < 114 GeV [SB]- So, the LEP limit constrains the low mass region for which there is 
no bound from the current direct detection. Furthermore, there is a strong limit from the 
2011 year data at the LHC on the Higgs-like scalar below 200 GeV, restricting the mixing 
angle to be as small as sin 2 6 = 0.1 [37]. Thus, close to the cancellation zone, the LEP and 
LHC bounds on the mixing angle are much stronger than the direct detection bound. On 
the other hand, away from the cancellation zone, the latter can be as strong as the former 
above m 2 = 300 GeV. 

On the other hand, for m 2 > 2m i or m 2 > 2M X , the singlet-like scalar can also 
decay into a pair of Higgs-like scalar or fermion dark matter with decay rates shown in 
appendix A. For a sizable mixing determined by the Higgs portal coupling Xhs and/or a 
sizable coupling between dark matter and the singlet scalar, the additional decay modes 
can easily dominate the SM Higgs decay modes of the singlet-like scalar: Xhs ^ 0.1 or 
\ x > 0.1 above the WW threshold [24"ll32"] . In particular, the dark matter coupling A x to 
the mediator can be sizable unless heavy fermion mass running in loops is not too heavy. 
In this case, the current LHC limit would not be directly applied. 



5 A UV complete model 

The toy model given in Section 2 contains dimension-5 interaction terms between the 
axion-like scalar and the electroweak gauge bosons. In this section, we consider the PQ 
symmetric extension of the MSSM0 as a UV completion of the toy model for generating the 
electroweak PQ anomalies. In this model, we can generate the electroweak PQ anomalies 
by the superparticles of the Higgs doublets, the Higgsinos, without introducing extra SM 

6 Different realizations of the NMSSM with PQ symmetry have been proposed recently [38] where SUSY 
breaking sector respects PQ symmetry. But, in our case, SUSY breaking sector breaks PQ symmetry while 
it respects discrete R symmetry. 
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non-singlets. In this UV completion, there are the invisible axion multiplet A and extra 
singlets, S and Xi,2, as well as two more singlets, X and Y, for the hidden sector SUSY 
breaking. 

Assuming that the singlets have nonzero PQ charges and R-charges as in Table 1, the 
superpotential for the extended Higgs sector is given by the following, 

W vis = X h SH u H d + X x S X iX2 • (22) 

The first term provides the fi term after the singlet S gets a VEV while the second term 
leads to the axion coupling to a Dirac fermion dark matter composed of Xi an d X2 for a 
nonzero singlet VEV. In our model, the QCD anomaly is generated only by heavy quarks 
coupled to the invisible axion multiplet A. PQ symmetry transformations are defined as 
A — > A + i8pQ and — > e tqi PQ ^fi for all the other chiral superfields with PQ charge q^ in 
the model. The model is similar to KSVZ axion models, because there is no direct coupling 
of the invisible axion to the Higgs doublets. 





Q,L 


17, iV 


D,E 


H u 


H d 


s 


Xi 


X2 


X 


Y 


PQ 


Qi 





-?i + <?2 


-Qi 




?1 + <?2 


Qx 


-Qx -Q1-Q2 


Qx 


2(gi + &) 




1 


1 


1 








2 








2 






Table 1: PQ and Zf charges. Right-handed neutrinos N are also included for neutrino 
masses. 



Discrete R symmetries provide a solution to the \i problem and guarantees the proton 
stability in the MSSM [39HU] . In particular, for S77(5)-GUT compatible R-charges, and 
Zg are the only discrete R symmetries that allow for a singlet extension of the MSSM [£D] . 
A cubic interaction and/or linear term for the singlet S is consistent with those discrete 
R symmetries but they are forbidden by PQ symmetry in our model. A non-perturbative 
dynamics in hidden sector might generate the \x term and other PQ breaking singlet mass 
terms in the superpotential. But, in our model, we assume that the dominant PQ symmetry 
breaking in the Higgs sector comes from a spontaneous PQ breaking at high scale. 

After the singlet S gets a VEV, the singlet fermions, xi an d X2, obtain a Dirac mass 
as M x = X x (s). If this singlet Dirac fermion is the lightest superparticle, it can be dark 
matter with odd R-parity as in the toy model in Section 2. Then, the Higgsinos generate 
the electroweak anomaly couplings to the axion coming from S as follows, 

ci = c 2 = ~(qi + 12)- (23) 

Now we discuss a microscopic model for SUSY breaking to generate the soft mass 
parameters. Assuming that SUSY breaking multiplets X, Y have PQ charges qx and 
qy = 2(<7i + q 2 ), respectively, the effective superpotential for the SUSY breaking sector is 

W hid = filXe- qxA + W {1 + li\Ye-^ A ). (24) 
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Here, we have introduced the R-symmetry breaking in terms of a constant superpotential 
Wo- After the saxion is stabilized, the above superpotential leads to nonzero F-terms, 
Fx = n\e~ qxA and Fy = W^[i\e~ qYA . Therefore, SUSY breaking sector also breaks PQ 
symmetry after the scalar partner of the axion multiplet is stabilized [2]. Henceforth we 
take two F-terms to be comparable. Furthermore, we can write the PQ and R invariant 
effective interactions composed of the invisible axion and SUSY breaking fields as follows, 

/ d4d ^d Y ^ S2 + J d2 ° aXe ~ (qx ~ qi ~ q2)AH ^ H d + f d A 6^Y^ Xe- {qx ~ qi - q2)A S + h.c. (25) 

with M being the messenger scale. The first term corresponds to a supersymmetric singlet 
mass term of order scalar soft mass from nonzero Fy term as follows, 

1 F f 
AW vis = -fi s S 2 , Vs = ajt. (26) 

This supersymmetric singlet mass is crucial to make the extra singlet fermion to be heavier 
than the Dirac fermion dark matter. On the other hand, from eq. (1251) . the second term 
generates a B-term for Higgs doublets while the third term generates a singlet linear soft 
mass term. We note that the soft trilinear term is also generated by SUSY breaking with 
f d 2 9jjY e~ qvA SH u H d but they are suppressed as compared to soft scalar masses of order 
^~m~ ~ a*' because the SUSY breaking fields carry nonzero PQ charges. Nonetheless, 
gravity mediation can generate soft mass terms of order gravitino mass corresponding to 
all the terms present allowed in the superpotential. Therefore, we get the soft SUSY 
breaking terms as follows, 

-Aoft = mljH u \ 2 + mljH d \ 2 + ml\S\ 2 + ml i \xi\ 2 + m 2 X2 \ X 2\ 2 + {^B sf i S S 2 + h.c. 

+ {X h A x SH u H d + X x A x S X iX2 + B^H u H d + B xf i X iX2 + B s m 2 S + h.c.) (27) 

where A\ is of order gravitino mass and 

~ B x fi ~ aF x e- {qx ~ qi ~ q2)A , B s m 2 ~ e~ {qx - qi ~ q2)A (28) 

If the above B-terms are of order the soft scalar mass, for \Fx\ ~ |^V|> we need 



7 A possible domain- wall problem could arise after R-symmetry breaking of order 10 11 GeV related to 
Wo in gravity mediation, but it depends on the reheating temperature after inflation. As far as the bound 
on the reheating temperature coming from gravitino problem, about 10 9 GeV, is satisfied, there is no 
domain wall produced after inflation, along the line of discussion in Ref. [39ll40j . 
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PQ breaking linear soft mass for the singlet stabilizes the singlet at a nonzero VEV 
and gives the S axion a nonzero mass. The minimization of the potential leads to [12] 

• O o (2A X + X hf i s )(s) 

sm2 ^ = " 0S= " <s> ' (30) 



^ ~~ -.2 I D , ,,2 I \2„2 ■ ^ 



Afatyu^AA + 2/xg) - B s m 2 
m 2 s + B s^s + (*s + ^ 2 h v 2 
The scalar potential in NMSSM is more predictive because the Higgs quartic coupling is 
given by the gauge coupling and there is no singlet self-coupling. However, as in the toy 
model, a mixing between CP-even Higgs and singlet is possible due to the singlet coupling 
to the Higgs doublets. 

The A-term in eq. (127)1 and the F-term for the S singlet from the superpotential, 
W vis + AW vis , can mix between the pseudo-scalar Higgs and the S axion, so the tree-level 
DM annihilation into a pair of SM particles through the mixing could be large. Thus, in 
order to explain the gamma-ray line with the correct relic density, the mixing between the 
pseudo-scalar Higgs and the S axion should be suppressed. From the gamma-ray constraint 
that the extra tree-level axion mediation through the mixing is smaller than the one-loop 
induced counterpart, if the Higgs pseudo-scalar and S axion masses, i.e. uia and m a , are 
comparable, we need to impose the soft mass parameter as 



1^ -"*!<, _^_„,/ ; ^L K o.01. (32) 



m a ~ V 327r 2 t> s (m a + M x ) V 9Qtt 2 v s 

If the pseudo-scalar Higgs mass is much larger than the S axion mass, the amount of a 
tuning on \A\ — fis\ is reduced by a factor of ^ SL . Then, the S axion couples dominantly 
to the electroweak gauge bosons in the SM through the anomalies, playing a role of the 
mediator between dark matter and the SM. Ignoring the mixing of pseudo-scalars for 
A\ ~ /Us, we obtain the pseudo-scalar masses as 



2 2(B^ + A x (s) + \ h fx s ) foo\ 

ml = j-r (\(A X + ti s )v u v d - Bsm 2 } - 2B S fi S - (34) 

We note that there are extra fields from the dark matter and messenger sectors in 
the supersymmetric models as compared to the toy model in Section 2: an extra Higgs 
doublet, the scalar partners of Dirac fermion dark matter and the fermionic partner of 
the S singlet. But, the scalar partners of Dirac dark matter are not relevant for PQ 
and electroweak symmetry breaking. Furthermore, the extra singlet superparticles can be 
heavier than dark matter with mass M x so the additional annihilation channels into a 
singlet pair can be kinematically suppressed. Moreover, since the extra CP-even Higgs is 
heavier than the SM-like Higgs, the DM annihilation cross section is determined dominantly 
by the mediation channels with the S axion and lighter CP-even scalars as discussed in 
our toy model in Section 2. 
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6 Conclusion 



We have considered a Dirac singlet fermion as dark matter that communicates with the 
SM by a complex scalar mediator. Identifying a U(l) global symmetry in the dark sector 
with PQ symmetry, a spontaneous breakdown of PQ symmetry at high scale generates a 
soft PQ-breaking mass of weak scale for the axion part of the complex scalar in a CP- 
invariant fashion. After the complex scalar gets a VEV, the effective axion interactions to 
the electroweak gauge bosons are generated by anomalies in the presence of extra heavy 
fermions with axion coupling. This opens up a possibility that the axion mediates dark 
matter interactions close to the resonance, m a ~ 2M X , such that dark matter annihilates 
into 77 and Z'-f with sizable branching fractions while reproducing the relic density. If the 
Fermi gamma ray line is confirmed, there will be interesting signatures to be pursued for 
DM direct detection and Higgs-like scalar searches below 300 GeV at the LHC in a near 
future. 

We also have presented a ultra-violet complete model that accommodates the axion 
coupling to the electroweak gauge bosons naturally by the singlet coupling to Higgsinos in 
the NMSSM with PQ symmetry. The SUSY extension relies on the specific PQ breaking 
soft mass terms in the NMSSM that are obtained in the presence of a discrete R symmetry. 
It would be worthwhile to investigate the implications of the gamma ray line on the Higgs 
boson and SUSY searches in this context. 
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Appendix A: Decay rates 

In the presence of the anomaly interactions between the axion and the electroweak gauge 
bosons, cv 1 V2 ae fj.upaFy^ Fy^ , the decay rate for a — > V1V2 with two gauge bosons, V\ and 
V2, having masses M\ and M 2 , respectively, is 

T(a ->• V1V2) = — s v \c Vl v 2 I 1 2 ) i 1 2 ) ( A - X ) 

with sy being the symmetry factor for the final states, for instance, sy = N\ for N identical 
final states. Furthermore, if the axion mass is larger than twice the dark matter mass, the 
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axion can decay into a dark matter pair. Then, the total decay rate of the axion is given 
by 

r a = r a ( 77 ) + r a (z 7 ) + r a (zz) + r a (ww) + r a ( X x) (A.2) 

where 

777 

T a (77) = -f\c~ r( \ 2 , (A.3) 



M|\3 
27r 1 " " V m 2 



r a (^7) = U-^f , (A.4) 



7/1 



7T 



4M|\3/2 



r a (ZZ) = -±\ Czz \*[l-—£\ , (A.5) 



r.W - f\c WW Hl~^f\ (A.6) 
|AJ 2 / 4m 2 \ 1/2 

16?T V 772^ / 

Here, the anomaly couplings are related to the original parameters in eq. ([1]) as 

(ciai cos 2 W + c 2 a 2 sin 2 (A-8) 

(c 2 a 2 - cia l )sm(26 w ), (A.9) 

(c 2 a 2 cos 2 ^vk + ci«i sin 2 ^ly), (A. 10) 

(A.ll) 

We also consider the decay rates of CP-even scalars. They can decay into the SM 
particles due to the mixing with Higgs boson as the SM Higgs does. Each partial decay 
width of the CP-even scalars into an SM particle pair is obtained from the one of the 
SM Higgs multiplied by sin 2 9 for s and cos 2 9 for h. If kinematically allowed, the heavier 
CP-even scalar s can decay into a pair of the lighter CP-even scalar h and the CP-even 
scalars can decay into a dark matter pair by the direct coupling. So, the decay rates for 
the possible additional decay modes are 





IQttVs 




1 


cz 7 = 


16tiv s 




1 


czz = 


IQttVs 






c 2 a 2 


C\VW = 


87TV S 



l(s^hh) = ^LJl-^ (A.12) 

2\ 3/2 



nS^xx) = %j^cos 2 ^l-^) , (A.13) 



16?T \ 7772 

IX l 2 m, / Am 2 \ 3 ^ 2 

T(h^XX) = ^fpisin 2 ^ (l--^) • (A.14) 
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Appendix B: Dark matter annihilation cross section 

First, due to the dark matter coupling to the axion and the axion anomaly interactions, 
the cross section times relative velocity for XX V1V2 is given by 



32tt 



l\ |2| |2 
s V\^x\ \ c ViV 2 \ - 



m 



i _ (M, + M 2 ) 2 ^ 3/2 ^ ^ (M l -M 2 f ^ 3/2 

(B.l) 

with s being the center of momentum squared. Then, the velocity averaged cross section 
for dark matter annihilation with axion mediation is 



(w) a = (crv) yy + (crv)z-y + (crv) zz + (av) 



WW 



where 



(av)yy 

(<7V) Z1 

(crv)zz 
(av) ww 



1 


\K 2 


c 77 | 2 






167T 




m 2 a ) 2 + m 2 a T 2 a ' 


1 


|A X | 2 | 








32vr 




- m 2 a ) 2 + m 2 a T 2 a 


1 


\K 2 


c Z z\ 2 




A. 


167T 




~ ™ 2 ) 2 + ml^l 


1 


|A X | 2 | 


Cww\ 




16M* 

A, 


32vT 




- m 2 a ) 2 + m 2 a Y\ 



V AM 2 ) ' 



M 2 \ 3 / 2 



M 2 X ) 
- M 2 X 



(B.2) 

(B.3) 
(B.4) 
(B.5) 
(B.6) 



Second, dark matter can also annihilate into an SM particle pair by the Higgs portal 
interaction to the real scalar partner of the axion. The DM annihilations through the 
CP-even scalars are p-wave suppressed unlike the counterpart of axion mediation as shown 
below. The dominant annihilation cross section coming from the CP-even scalar interaction 
is composed of 

(av) s = (crv)f f + (crv) ww + (av)zz + (B.7) 
The partial annihilation cross section into an SM fermion pair is 



(av) 



ff 



3iV c |A x | 2 m 2 sin 2 (20) 



32nv 2 



P^ 2 M 2 ((m\ - m 2 2 ) 2 + (m 1 T 1 - m 2 Y 2 f) 



x 1 



m 2 \3/2 t 
M 2 ) M~ x 



(B.8) 



with T being the temperature of the universe and Pi = [(AM 2 — m 2 ) 2 + m 2 r|] 1 (i=l, 2). 
From the gauge-Higgs interactions, cyhV^V^ 1 with 



cw = —2 



M w 



cz 



M 2 



(B.9) 



18 



the annihilation cross section into a gauge boson pair is similarly obtained as follows, 



(av) vv = s v 



3|Axl 3 C 2 ^4 2(2g) (K - mif + (m iri - m 2 T 2 f) 

/ My 3 My \ / i^y/ar 



with sy being the symmetry factor. Finally, from the interactions between CP-even scalars, 
J,eii h 3 + s /i 2 , the annihilation cross section into a Higgs-like pair is 



\2 

(av) 



|A X , 



'hh 



16tt 



/ m 2 \ 1/2 T / \ 



where 



• 4 " = ' (R12) 

sin 2 6 M x (5M 2 - 2m 2 ) ^ . , 

A ts = x l * l) Y, Pi Re[A x a*] (4M 2 - m 2 ) , (B.13) 

V2(2M 2 -m 2 ) i=1>2 

A " = I ( 2 Pi|Si|2 + 2P i P 2Maia* 2 ] {(4M 2 - m 2 )(4M 2 - m 2 ) + m^r^} 

\i=l,2 , 

(B.14) 

with 

3m 2 

ai = (v s cos 3 + v sin 3 0) sin , (B.15) 

Vol? 



sin (20) (2m 2 + m 2 ) . „. .„ . 

a 2 = 1 — ^ 1 ^(vsin0-<u s cos0)cos0. B.16 

2t> ,f 
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